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To easily calculate statistical properties of pairs correlated through Schmidt decomposition, as 
commonly used in Quantum Information, we propose a "commutator formalism" for these single- 
index pairs, somewhat simpler than the one we developed for double-index Wannier excitons. We 
use it here to get the pair number threshold for bosonic behavior of A'' pairs through the requirement 
that their number operator mean value must stay close to A''. While the main term of this mean value 
is controlled by the second moment of the Schmidt distribution, so that to increase this threshold, 
we must increase the Schmidt number, higher momenta appearing at higher orders lead to choosing 
a distribution as flat as possible. 



Over the last decade, we have extensively studied [1] 
fermion pairs making semiconductor excitons. We have 
shown that excitons mainly interact through the Pauli ex- 
clusion principle between their fermionic components [2] . 
This Pauli blocking gives rise to a set of scatterings which 
correspond to carrier exchanges in the absence of carrier 
interaction. Being by construction dimensionlcss, these 
"Pauli scatterings" control nonlinear effects induced by 
unabsorbed photons, as easily understood from a bare di- 
mensional argument: scatterings involving Coulomb in- 
teraction are energy-like quantities; they must appear 
with energy denominators which for non-linear optical 
effects, are of the order of the photon detuning. Detun- 
ing of non-absorbed photons being by construction large, 
these Coulomb processes thus are negligcablc in front of 
pure fermion exchanges. In the same way, just because 
they arc dimensionlcss, these Pauli scatterings cannot ap- 
pear in effective hamiltonia ns for bosonized excitons [3] . 
As a result, such exciton effective hamiltonians miss a 
large amount of physical effects, whatever their effective 
exciton-exciton scatterings. 

Usual semiconductor excitons are Wannier excitons. 
They are made of free electrons and free holes. These 
double-index excitons i = {Qi,Vi) where Qj is the center- 
of-mass momentum and Vi the relative motion index, arc 
linear combination of double-index fermion pairs. Their 
creation operators read through their wave function in 
momentum space as 



(1) 



We have proposed a "commutator formalism" [1, 2] to 
easily handle the consequences of the Pauli exclusion 
principle between N Wannier excitons. Pauli scatterings 
^imi) fermion exchanges between excitons starting 
in states ans ending in states {m,n), in the ab- 

sence of Coulomb process, formally appear through two 



commutators, 

e{aC,) + a( 



= Smi - Drni and 

These Pauli scatterings arc the 

keys to explain and better predict physical effects involv- 
ing unabsorbed photons. They however are rather com- 



ni 



plcx quantities. 

Later on, we turned to Frenkel excitons [4]. Being 
made of ion-site excitations, these are single-index exci- 
tons. Their creation operators read 



1 



(2) 



where Ng is the number of ion sites located in a periodic 
lattice R„. The fact that the Frenkel exciton wave func- 
tion is just a phase, induces important simplifications on 
the consequences of Pauli blocking on these correlated 
pairs. 

More recently, wc turned to Cooper pairs which also 
arc linear combination of single-index pairs 



(3) 



However, problems raised in BCS superconductivity re- 
quire not only an exact treatment of the Pauli exclusion 
principle but also of the potential between up and down 
spin electrons, to possibly generate the singular potential 
dependence of the energy [5]. This is why the "commu- 
tator formalism" we have developped for this problem [6] 
stayed at the free pair level. 

There is another field in which Pauli blocking between 
composite bosons plays an important role: Quantum In- 
formation. In this field, the correlated pairs are usually 
written [7, 8] through their Schmidt [9, 10] decomposition 



B 



t _ 



(4) 



These single-index pairs thus have similarities with 
Frenkel excitons although the Schmidt distribution ^/X^ 
may not be flat, i.e., it can differ from just a phase. 
In view of the importance of the field and the difficulty 
to properly handle the Pauli exclusion principle through 
a brut force counting of the amount of blocked events, 
it is highly desirable to develop a "commutator formal- 
ism" appropriate to fermion pairs correlated through a 
Schmidt decomposition which are the relevant pairs in 
Quantum Information. 
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We wish to stress that, according to Schmidt theorem, 
any two-fermion state can be written through a Schmidt 
decomposition. The double-index Wannier exciton dc- 



fined in Eq.(l) would then read B] = y '•p,t^p,ir^p,i- 
However, the operators appearing in this decomposition 
depend on the state i of the exciton at hand. Since, in 
physically relevant problems dealing with Wannier ex- 
citons, the excitons are scattered into different states, 
the Schmidt decomposition of Wannier exciton creation 
operators, with different elementary fermion operators 
for each i exciton, is not appropriate to approach the 
Wannier exciton physics. This is why we have been led 
to develop a formalism in which the double-indices of 
these composite bosons are explicitly kept. Fermion ex- 
changes between excitons appear through Pauli scatter- 
ings A(^j), these dimensionless scatterings allowing Wan- 
nier excitons to change states from to {in,n). The 
introduction of Pauli scatterings however is unnecessary 
for single-index pairs. The formalism we here propose 
is far simpler. It is based on a convenient mathemati- 
cal quantity, the "generalized correlated pair" defined in 
Eq.(5). Using it, it is possible to rederive within a few 
lines, some important Pauli blocking results obtained in 
the past through far heavier procedures. 

In this Letter, we first develop a formalism which al- 
lows an easy handling of Pauli blocking between N single- 
index correlated pairs. We then use it to calculate the 
mean value of the pair number operator as well as the 
variance of this number. In a last part, we discuss the 
appropriate shape of the Schmidt distribution to have a 
pair number mean value as close to N as possible. 

Formalism: We consider free fermion pairs charac- 
terized by a single index /. Their creation operators 



read Bj 



where a 



and bj 



are creation opera- 



tors for their fermionic components. The s/'s are such 
that [Bl,Bj] = while [Bi,,bI] = 5vi - where 

Di'i = Si'i{ajai + blbi). From them, we construct a set of 
"generalized correlated pair" operators as 



(5) 



with n = (0, 1, 2, . . . ) and (pi normalized by J^i I'Pfl = 1- 
These operators are such that [C,^„, C,t ] = 0, while 



(6) 



Tra is a scalar equal to the (m + l)-moment of the ipi 
distribution, namely 



while A„ = Ei \<t>'i\"'^Ha\ai + b\bi). Since [aJa/,B/,] ^ 
Sii'Bj, = [blbi,Bl], we readily find 

[Dm,C}j] — 2Cj^^^_^_i (8) 



For N correlated pairs with creation operators Cq, it 
is convenient to note that the iteration of Eq.(8) leads to 

[Dm,CQ ] = [-Dm, CqJCq +CQ[Dfn,CQ ] 

= 2iVC,U,Cr-^ (9) 

So, using Eq.(6), we get by iteration 

r'^N] — \r< /^ti/^t-'v-i I ^tf/^ ^t-'v-ii 

= 7VCr-'(T-™ - An) - N{N l)Cr-'cl^, (10) 

As evidenced below, all statistical properties of N cor- 
related pairs with creation operator Cq follow from 
this commutator. For the (unnormalized) iV-pair state 
\iPn) ~ Cq^|0), it in particular gives, since Dm\Q) = 

Cmli'N) = NT„-^^PN-l) - N{N - l)Cl^^\^PN~2)■ (H) 

Number mean value for N correlated pairs: Let 

us write the norm of the \iPn) state as 

{MM = (0|Co^C(l^|0) = iVIFjv (12) 

_F/v is a crucial quantity [11, 12] for many-body effects in- 
duced by the non-bosonic behavior of composite bosons, 
"cobosons" in short. Fjq would be equal to 1 for C\ cre- 
ating an elementary boson. For cobosons, i^^r, still equal 
to 1 for N = (0, 1), decreases when N increases because 
more and more pair states are excluded from the Cq sum 
due to Pauli blocking: this is the so-called "moth-eaten 
effect". 

We want to determine the mean value of the pair num- 
ber operator N = CqCq in the |V'7v) state. To do so, it 
is convenient to first note that, as CgC|=C|CQ while 
To = 1, Eq.(ll) used for N and N+l gives [13] 



ClCo\'<i^N) = N\^n) - NiN ~ l)Cl\M~i) 
N - 1 



1^ 



iV+ 1 



Co IV. 



'JV+l/- 



(13) 



From it, we readily get the pair number mean value [11, 
14] in the A'^-pair state [f/'Af) as 

(iV)^^iM^^l+(^_l)^|±i=Ar(i_c^.) (14) 



{■)pN\tp. 



N 



Fn 



where is the fraction of composite bosons which devi- 
ates from elementary bosons: The smaller C^v, the closer 
to an elementary boson behavior, ^jv precisely reads 



c 



N 



N - 1 
N 



(1 



N + l . 



F 



N 



(15) 



This "deviation fraction" reduces to zero for elementary 
bosons, i.e., Fn = 1 whatever TV. It also reduces to 
zero for iV = 1: two composite particles are needed to 
evidence their statistical nature. Note that ^ Civ ^ 1 
since, due to the "moth-eaten effect", F^ is a decreasing 
function of N. 
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Using Eq.(13), we also find the variance of tlic particle 
number as 



{N)n 

(iv-i)2 n- 



iV+l 



1+(A^-1) 



Fn + 1 



Fn+1 
Fn 



(16) 



We check that, as (n, this variance reduces to for ele- 
mentary bosons, Fn = 1, and also for iV = 1. 

When N increases, we physically expect cobosons to 
" shrink" more and more due to the Pauli exclusion prin- 
ciple between their fermionic components: cobosons be- 
come more and more difTerent from a set of single co- 
bosons. As a result, the fraction of cobosons which devi- 
ates from an elementary boson behavior should increase 
with N from Ci = to Ca? ~ 1 above a certain pair 
number threshold. In order to determine the N scale 
over which this occurs, we must explicitly calculate the 
Fn+iI Fn ratio which indeed is the key parameter in this 
problem. 

Normalization factor ratio. By writing (V'atI as 
(V-AT-iICo , Eq.(ll) gives, since tq^I 



('^ArlV'Ar) = N {ipN -l\i>N -l) 

This equation, used for - 



-N{N-l){^N-l\C\\4^N-2)- 

(17) 

1, yields 



F; 



N+1 



1 - N 



(18) 



Fn {iPn\iPn) 

This already shows that when the distribution is flat, 
= exp{iipi)/y/l^ for 1 < I < N* , with N* being the 

"Schmidt number" and 0; = otherwise, Cj reduces to 

Cl/N*, so that 



r Fn+1 \ 
I Fn > fia 



N 

flat ^ N* 



(19) 



As a result, Fn+i = for A^ ^ A^*: the state \iI.'n+i) 
then reduces to zero. 

In the case of an arbitrary distribution, we can cal- 
culate the Cj matrix element in Eq.(18) using Eq.(ll). 
More generally, this equation gives 



{lljN-l\Cm\tpN)=NTjn{lpN-l\lpN-l) 

-N{N~l){^N-l\CUi\i^N^2) 



(20) 



Using Eq.(17) for N{t(;N-i\il'N-i), the above equation 
allows us to rewrite the Fn ratio as 



F, 



N+l 

'f^ 



+N^{N - 1) 



{lpN-2\C2 - TiCi\lpN-l) 



{tpN\^. 



'N 



(21) 



For a flat distribution, Cl^ = Cl/N*"' and r,„ = 1/A^*"; 
so, the last term in the RHS cancels, in agreement with 
Eq.(19). 



To go further, we iterate the process using Eq.(19). 
This gives 



Fn 



+1 



N 



1 - Nti + N{N - 1){t2 - rf) 



Fn-'. 
Fn 



Ar^(A-l)^(A-2)i^:^^^i^f:Il4^^ (22) 



(V'AflV'Af) 

And so on ... Iteration using Eq.(20) leads to 



"N+l 



F 



N 



l-ATl+Ar(A^-l)|(T2-T^) 

-(A^-2)(T3-T2ri; 



F, 



N-2 

Fn 

Fn-3 

Fn 



+ (Af-2)(iV-3)(r4-r3Ti) 



Fn-4 



N 



(23) 



For a fiat distribution, r„ = (l/A^*)"; the prefactors of 
the Fn ratios in the RHS cancel: and we recover Eq.(19). 
For a general distribution, these prefactors have alternate 
signs. Indeed, since tq = 1 

I I' I I' 

= ^Ei'^'ii'^?4''^''^''^'''}{''^'"'~''^''"'}^° ^^^^ 

This shows that the flatter the distribution, the smaller 
the corrections to the main term (1 — Nti). 

In view of the value for a flat distribution, dimen- 
sional arguments lead to scaling as (1/A"e//)™ where 
1/ri is the "effective Schmidt number" of the dis- 
tribution at hand. For A^ small compared to this number, 
Eq.(23) gives the Fn ratio as [15] 



F, 



N+l 



F, 



N 



1-Nti + N\t2-t^)-N'^{t3-3t2Ti+2tI)-- 



(25) 

For arbitrary A^, it is possible to show that the sum of 
terms in N{N—1) in Eq.(23) gives a positive contribution 
to the Fn ratio, so that this ratio is larger than (1 — 
Nti) in agreement with the beautiful inequality recently 
derived by Wootters's group [16], namely 



1 AT ^ Fn+1 , -, 

1 - Ntx ^ — — < 1 - n 

tN 



(26) 



To show it, we come back to Eq.(21). By noting that 
the C\ matrix element in Eq.(17) must be real, we can 
rewrite the A^(A^ — 1) factor in Eq.(21) as 

iV(V'jv-2|C2|V^w-i) - Afri(^jv-i|CiVjv-2) (27) 

We then replace Nti{iPn-i\ according to Eq.(9). This 
leads to 

l<>i^N-2\C2^\-^N-l) - A(A - \)CUjN-2) 

~l,^N\Cf\-^N-2) (28) 
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in which we again use Eq.(ll) to replace the bracket by 
ColV'Af)- As a result, Eq.(21) also reads 



Fn+1 
F, 



= 1 - TVn 



N 



The simplest way to show that the last term is positive 
is to expand it on free pair operators. We then find that 
the distribution appears through 



(#) (#) 

, (#) f ^2 

= 2 E l-^'l-'-lCl l-^fJ-l^l ^0 (30) 

the sums being taken over differcnts {h ■ ■ ■ In)- As a re- 
sult, Fn+i/Fn is larger or equal to f — Nti. 

To derive the upper bound in the inequality (26), we 
can also use Eq.(29). We then have to show that 

{N ~1)ti{iPn\^m) 

+N{N -1)\{^n\C\^\^n-2) - {'4'n~i\C2\^n)\ (31) 



is positive. Using Eq.(20) for ti{%1]m\'4'n) and noting that 
the 6*2 matrix element is real as seen from Eqs. (18,21), 
the above quantity also reads 



N -I 
N + 1 



{^N\Ci\iJN+i) + N{N - l)(Vjv|Cf IVjv-2) (32) 



This gives the first Fat's as Fi=ro=l, F2=l— ti, 
^3=1— 3ti-|-2t2, the general expression of Fjv reading a 



Fn = 1-A^(7V-1) ^ll-{N^2)^l2+{N-2){N~3)^l3 ■ 



which is the sum of two positive terms as easily seen by 
expending them on free pair operators. We wish to note 
that the upper bound in the inequality (26) just corre- 
sponds to Fn^i/Fn being a decreasing function of N 
since F2/F1 = (1 — n) as seen from Eq.(21). Actually, 
this TV independent upper bound leads to a "deviation 
fraction" larger than ti(1 — 1/N) which is not re- 
ally useful to determine the N threshold above which 
cobosons substantially deviate from bosonic behavior. 

Although the above formalism provides a very direct 
way to reach the Fn^i/Fn ratio relevant in the problem 
we here consider, we wish to mention that it is possible 
to recover this ratio through the general equation which 
links the various i<jv's. A simple way to get this equation 
which does not use the Pauli scatterings A(^j) introduced 
in previous derivations, is to start with Eq.(17) and use 
Eq.(20). We then get 

{^n\^n) = N{4,N-l\4'N-l) - N{N - 1)''ti{iPn-2\^N-2) 

+NiN ~ 1)\N - 2)2(V'^_3|C2|Vjv-2). (33) 

Iteration leads to an equation between the Fpf^s only 
which reads as 



Fn - FN-i-{N-l)TiFN-2+iN-l){N^2)T2FN^3- 



(35) 

with fii = ri/2, fi2 = T2/3, fj,3 = T3/4: - rf/S, and so 
on... as can be checked by inserting Eq.(35) into Eq.(34). 
This shows that F^ has an overextensive dependence on 
N, with terms in Nf{N/N*j:j), this overextensive depen- 
dence however disappearing in physical quantities which 
depends on Fn ratios only. 

"Deviation fraction" We now come back to the "de- 
viation fraction" defined in Eq.(15) and consider a few 
particular distributions. 

(i) Flat distribution: Using the value of the F^+i/Fn 
ratio for a fiat distribution given in Eq.(19), we readily 
find that the "deviation fraction" then takes a compact 
form, namely 



(36) 



We see that this "deviation fraction" Cjv increases when 
N increases from its bosonic value Ci = to C,m' = 
(1 - 1/A^*) and then stays equal to 1 for iV > {N* + 1), 
all cobosons differing from elementary bosons above this 
threshold. Eq.(36) shows that, for a flat distribution, the 
larger the Schmidt number N* , the larger the number 
of cobosons which possibly behave as elementary bosons. 
If we now turn to the variance of the particle number 
defined in Eq. (16), we find that, for a flat distribution, 
it reduces to zero for all N ^ as in the case of elementary 
bosons 



flat 



0. 



(37) 



(ii) Peak or canyon distribution: We now consider a 
distribution with the same Schmidt number N* , i.e., 
7^ for 1 ^ / ^ A^* but with |(/);| now peaked or 
depressed on I ~ Ic- It is clear that in this case too, 
Pauh blocking imposes C(J^|0) = for A^ > A^*. The 
question is to know if, for a given A^*, a peak distribution 
would help to decrease the fraction of A^ cobosons which 
deviates from bosons. Intuitively, the answer should be 
no, because a peak distribution can be approximated by 
a step function which fundamentally corresponds to a 
reduction of the Schmidt number. 

To illustrate it, let us consider a very simple case: all 
the 1 0ip's are equal except one, namely 10; p = {l—x)/N* 
for all I Ic in the {1,N*) range, tq = 1 then imposes 
= X -I- (1 — x)/N*; so, positive x's correspond to a 
peak and negative x's to a canyon. The fact that |(/>;^;^p 
and p must be both positive restricts x to — 1/(A^* — 
1) < X < 1. The X upper bound corresponds to only 
having the Ic state populated while the x lower bound 
corresponds to having this state empty. The second 
and third momenta of this distribution are given by 



(34) 



N*Tl 
N*^T2 



x^{N* - 
-3x^{N* 



I) 



l) + x^{N* ~ 1){N* - 2). (38) 
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Eqs. (15,38) then give the "deviation fraction" for two 
correlated pairs as 



C2 



Tl - T2 



1 



1 + (iV* - 2) 



1 



(39) 



We recover the flat distribution result of Eq.(36) for x = 
0, the above equation then giving C2 = 1/-/V*, and also 
for X = —1/{N* — 1) which corresponds to = 0: 

we then only have [N* — 1) states in the entanglement; 
so, C2 = 1/(A^* - 1). Eq.(39) shows that, in the range of 
physical x's, the fraction of two cobosons which deviates 
from elementary bosons is minimum for x — Q, i.e., when 
the \4>i\^ distribution is flat as physically guessed: to de- 
viate from a flat distribution, by either making a narrow 
hole (x < 0) or a narrow peak {x > 0) tends to increase 
this fraction. 

(Hi) Arbitrary distribution: We now turn to an arbi- 
trary distribution. The t„ momenta decrease from tq = 1 
when n increases. The broader the distribution, the 
faster the decrease, r„ scaling as (l/iV*j^)" where N*^j is 
the "effective Schmidt number" . To possibly have a large 
number of cobosons behaving as elementary bosons, the 
distribution must be broad enough for t„ to be much 
smaller than 1 for n = 1 already. For ^ 1 but still 
small enough for N/N*^^ to be small, Eq.(25) gives the 
"deviation fraction" as 



Cat « NTl-N^{T2-Tl)+N''{T^-'dT2T^+2Tl)^ 



(40) 



We see that, since t„ oc (1/A^*^j)", the first term of C,n 
scales as N/N*jj, the second term as {N/N*jjY ^-i^d so 
on... This shows that the scale for the number of co- 
bosons possibly behaving as elementary bosons is essen- 
tially related to the effective Schmidt number through 
the second moment ti of the Schmidt distribution, of- 
ten called purity: the larger -/V*^^ = 1/ti, the larger the 
correlated pair number for which C^r <c 1. This physi- 
cally leads to take the largest possible Schmidt number, 
corrections to the first term of Eq.(40), which depend on 
higher momenta, leading us to choose a 0; distributin as 
flat as possible. 

To conclude, we have constructed a "commutator 
formalism" for pairs correlated through the Schmidt de- 
composition which are the relevant pairs in Quantum In- 
formation. Through its key equation (11), this formal- 
ism allows us to recover results linked to Pauli block- 
ing obtained previously in various different contexts ( 
Wannier excitons [1], Frenkel excitons [4], Cooper pairs 
[6] and also Quantum Optics [13]) through far heavier 



procedures. This formalism is dcfinitly valuable for the 
reader not to be forced to go through more complicated 
derivations which only are of interest for problems dealing 
with double-index correlated pairs. The present formal- 
ism makes use of a " generalized correlated pair" creation 
operators Cjj which are convenient mathematical quan- 
tities for easy calculations cnvolving N correlated pairs 
with same creation operator Cq . 

We first use this " commutator formalism" to calculate 
the ratio of normalization factors F^^i/Fm in terms of 
the various momenta of the Schmidt distribution, this 
ratio being the relevant quantity for physical effects in- 
duced by the Pauli exclusion principle between N iden- 
tical correlated pairs. This formalism allows us to also 
rederive the beautiful inequality for the Fn^i/Fn ratio 
recently obtained by Wooters's group [16] through a care- 
ful counting of the amount of terms left after Pauli block- 
ing in combination of scalar products like (?/'7v|V'Jv)'s. 

Wc finally use this formalism to determine the Scmidt 
distribution for which N cobosons would be close to N 
elementary bosons, by looking at the mean value of their 
number operator. We find that for a flat distribution, i.e., 
for a distribution which is just a phase factor, the frac- 
tion of N cobosons which deviates from an elementary 
boson behavior takes a very compact form, in terms of 
the second moment of the Schmidt distribution, the vari- 
ance of this particle number then reducing to zero. In the 
general case, this "deviation fraction" has a more com- 
plex form in which enter higher momenta of the Schmidt 
distribution. The study of a "peak" or "canyon" distri- 
bution shows in a transparent way that, while an increase 
of the Schmidt number always is favorable, the flatter the 
distribution, the larger the amount of cobosons behaving 
as elementary bosons. 
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